Inhomogeneous Gilbert damping from impurities and electron-electron interactions 
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We present a unified theory of magnetic damping in itinerant electron ferromagnets at order 
including electron-electron interactions and disorder scattering. We show that the Gilbert damping 
coefficient can be expressed in terms of the spin conductivity, leading to a Matthiessen-type formula 
in which disorder and interaction contributions are additive. In a weak ferromagnet regime, electron- 
electron interactions lead to a strong enhancement of the Gilbert damping. 

PACS numbers: 76.50.+g,75.45.-t-j,75.30.Ds 



Introduction - In spite of much effort, a complete 
theoretical description of the damping of ferromagnetic 
spin waves in itinerant electron ferromagnets is not yet 
available.^ Recent measurements of the dispersion and 
damping of spin-wave excitations driven by a direct spin- 
polarized current prove that the theoretical picture is in- 
complete, particularly when it comes to calculating the 
lincwidth of these excitationsJ^ One of the most impor- 
tant parameters of the theory is the so-called Gilbert 
damping parameter which controls the damping rate 
and thermal noise and is often assumed to be indepen- 
dent of the wave vector of the excitations. This assump- 
tion is justified for excitations of very long wavelength 
(e.g., a homogeneous precession of the magnetization), 
where a can originate in a relatively weak spin-orbit (SO) 
interaction^. But it becomes dubious as the wave vector 
q of the excitations grows. Indeed, both electron-electron 
(e-e) and electron-impurity interactions can cause an in- 
homogeneous magnetization to decay into spin-flipped 
electron-hole pairs, giving rise to a contribution to the 
Gilbert damping. In practice, the presence of this contri- 
bution means that the Landau-Lifshitz-Gilbert equation 
contains a term proportional to — m x V^dtm (where m 
is the magnetization) and requires neither spin-orbit nor 
magnetic disorder scattering. By contrast, the homoge- 
neous damping term is of the form mxdtm and vanishes 
in the absence of SO or magnetic disorder scattering. 

The influence of disorder on the linewidth of spin 
waves in itinerant electron ferromagnets was discussed in 
Refs. HISS and the role of e-e interactions in spin-wave 
damping was studied in Refs. for spin-polarized liq- 
uid He^ and in Refs. [ToUTll for two- and three-dimensional 
electron liquids, respectively. In this paper, we present 
a unified semiphenomenological approach, which enables 
us to calculate on equal footing the contributions of dis- 
order and e-e interactions to the Gilbert damping pa- 
rameter to order q^. The main idea is to apply to the 
transverse spin fluctuations of a ferromagnet the method 
first introduced by Mermin^^ for treating the effect of 
disorder on the dynamics of charge density fluctuations 
in metals.-iS. Following this approach, we will show that 
the q^ contribution to the damping in itinerant electron 
ferromagnets can be expressed in terms of the transverse 



spin conductivity, which in turn separates into a sum of 
disorder and e-e terms. 

A major technical advantage of this approach is that 
the ladder vertex corrections to the transverse spin- 
conductivity vanish in the absence of SO interactions, 
making the diagrammatic calculation of this quantity a 
straightforward task. Thus we are able to provide explicit 
analytic expressions for the disorder and interaction con- 
tribution to the q^ Gilbert damping to the lowest order 
in the strength of the interactions. Our paper connects 
and unifies different approaches and gives a rather com- 
plete and simple theory of q^ damping. In particular, we 
find that for weak metallic ferromagnets the q^ damping 
can be strongly enhanced by e-e interactions, resulting in 
a value comparable to or larger than typical in the case 
of homogeneous damping. Therefore, we believe that the 
inclusion of a damping term proportional to q^ in the phe- 
nomenological Landau-Lifshitz equation of motion for the 
magnetization^^ is a potentially important modification 
of the theory in strongly inhomogeneous situations, such 
as current-driven nanomagnetsSi and the ferromagnetic 
domain-wall motioi*i^Ji 

Phenomenological approach - In Ref. Mermin con- 
structed the density-density response function of an elec- 
tron gas in the presence of impurities through the use 
of a local drift-diffusion equation, whereby the gradient 
of the external potential is cancelled, in equilibrium, by 
an opposite gradient of the local chemical potential. In 
diagrammatic language, the effect of the local chemical 
potential corresponds to the inclusion of the vertex cor- 
rection in the calculation of the density-density response 
function. Here, we use a similar approach to obtain the 
transverse spin susceptibility of an itinerant electron fer- 
romagnet, modeled as an electron gas whose equilibrium 
magnetization is along the z axis. 

Before proceeding we need to clarify a delicate point. 
The homogeneous electron gas is not spontaneously fer- 
romagnetic at the densities that are relevant for ordinary 
magnetic systems;^'^ In order to produce the desired equi- 
librium magnetization, we must therefore impose a static 
fictitious field Bq. Physically, Bo is the "exchange" field 
Bex plus any external/applied magnetic field Bq^^ which 
may be additionally present. Therefore, in order to calcu- 
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late the transverse spin susceptibility we must take into 
account the fact that the exchange field associated with 
a uniform magnetization is parallel to the magnetization 
and changes direction when the latter does. As a result, 
the actual susceptibility XQf,(q, cj) differs from the sus- 
ceptibility calculated at constant Bq, which we denote 
by Xabi'^,'-^), according to the well-known relationjii 



Xabi^^^)=^Xali^^'^)-^^ab- (1) 

Here, Mq is the equilibrium magnetization (assumed to 
point along the z axis) and uj^x — iBex (where 7 is the 
gyromagnetic ratio) is the precession frequency associ- 
ated with the exchange field. Sab is the Kronecker delta. 
The indices a and b denote directions {x or y) perpen- 
dicular to the equilibrium magnetization and q and uj 
are the wave vector and the frequency of the external 
perturbation. Here we focus solely on the calculation of 
the response function x because term uJexSab/Mo does 
not contribute to Gilbert damping. We do not include 
the effects of exchange and external fields on the orbital 
motion of the electrons. 

The generalized continuity equation for the Fourier 
component of the transverse spin density Ala in the di- 
rection a {x or y) at wave vector q and frequency lo is 

-iujMaiq, uj) =- i-fq ■ ja(q, uj) - uJoeabMbiq, lu) 

+ jMoeabB^^''iq,Lu), (2) 

where _Bj^PP(q, oj) is the transverse external magnetic field 
driving the magnetization and loq is the precessional fre- 
quency associated with a static magnetic field Bq (in- 
cluding exchange contribution) in the z direction, ja is 
the ath component of the transverse spin-current density 
tensor and we put h — 1 throughout. The transverse 
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Levi-Civita tensor eab has components e^x 
£xy = ~^yx = 1, and the summation over repeated in- 
dices is always implied. 

The transverse spin current is proportional to the gra- 
dient of the effective magnetic field, which plays the role 
analogous to the electrochemical potential, and the equa- 
tion that expresses this proportionality is the analogue of 
the drift-diffusion equation of the ordinary charge trans- 
port theory: 



ja(q,t^) = iqa±_ 



Ma{q,uj) 



(3) 



where a±^ (= Uxx or Uyy) is the transverse dc (i.e., lo — Q) 
spin-conductivity and x± = Mq/luq is the static trans- 
verse spin susceptibility in the g — > limit. Just as in 
the ordinary drift-diffusion theory, the first term on the 
right-hand side of Eq. ([3]) is a "drift current," and the 
second is a "diffusion current," with the two canceling 
out exactly in the static hmit (for q ^ 0), due to the 
relation Ma(0,0) = 7x^5^^(0,0). Combining Eqs. © 
and ^ gives the following equation for the transverse 



magnetization dynamics: 



-lUJOab H r dab + (^Q(^ab ] Alb 



{MoCab + lCf^q^6ab) lB\ 



which is most easily solved by transforming to the 
circularly-polarized components M± = Aix ± iM^, in 
which the Levi-Civita tensor becomes diagonal, with 
eigenvalues ±i. Solving in the "-f " channel, we get 



Mo — i^a_\_Q^ 
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from which we obtain to the leading order in uj and 
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(6) 



The higher-order terms in this expansion cannot be legit- 
imately retained within the accuracy of the present ap- 
proximation. We also disregard the q^ correction to the 
static susceptibility, since in making the Mermin ansatz 
([5]) we are omitting the equilibrium spin currents respon- 
sible for the latter. Eq. ([6]), however, is perfectly ade- 
quate for our purpose, since it allows us to identify the 
q^ contribution to the Gilbert damping: 
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^mx+-{q,uj) J<J±q'^ 
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(7) 



Therefore, the Gilbert damping can be calculated from 
the dc transverse spin conductivity (7j_, which in turn 
can be computed from the zero-frequency limit of the 
transverse spin-current — spin-current response function: 



1 



jj^^ Sm((^^^]^ SjaPia] Sia'Pia))uj 



where Sta is the x or y component of spin operator for the 
ith electron, pia is the corresponding component of the 
momentum operator, is the effective electron mass, V 
is the system volume, N is the total electron number, and 
{{A]B))^ represents the retarded hnear response func- 
tion for the expectation value of an observable A under 
the action of a field that couples linearly to an observable 
B. Both disorder and e-e interaction contributions can 
be systematically included in the calculation of the spin- 
current — spin-current response function. In the absence 
of spin-orbit and e-e interactions, the ladder vertex cor- 
rections to the conductivity are absent and calculation 
of o-j_ reduces to the calculation of a single bubble with 
Green's functions 



G'Ta(P'^) 
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ep + £_F±tJo/2 + i/2T' 



(9) 



where the scattering time in general depends on the 
spin band index s =T,i- In the Born approximation. 
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the scattering rate is proportional to the electron den- 
sity of states, and we can write | = tv jv^ \^^ where Vg 
is the spin-s density of states and v — (y\ + v\)l1. r 
parametrizes the strength of the disorder scattering. A 
standard calculation then leads to the following result: 



_dis 
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(10) 



This, inserted in Eq. ([7]), gives a Gilbert damping pa- 
rameter in full agreement with what we have also calcu- 
lated from a direct diagrammatic evaluation of the trans- 
verse spin susceptibility, i.e., spin-density — spin-density 
correlation function. From now on, we shall simplify the 
notation by introducing a transverse spin relaxation time 



(11) 



where Eps — m^,Vpg/2 is the Fermi energy for spin-s elec- 
trons and n is the total electron density. In this notation, 
the dc transverse spin-conductivity takes the form 
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(12) 
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Electron- electron interactions ~ One of the attractive fea- 
tures of the approach based on Eq. ^ is the ease with 
which e-e interactions can be included. In the weak cou- 
pling limit, the contributions of disorder and e-e inter- 
actions to the transverse spin conductivity are simply 
additive. We can see this by using twice the equation of 
motion for the spin-current — spin-current response func- 
tion. This leads to an expression for the transverse 
spin-conductivity ([8]) in terms of the low-frequency spin- 
force — spin-force response function: 
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.(13) 



Here, Fia — Pia is the time derivative of the momentum 
operator, i.e., the operator of the force on the ith elec- 
tron. The total force is the sum of electron-impurity and 
e-e interaction forces. Each of them, separately, gives a 
contribution of order jweip and |feeP, where v^i and Vee 
are matrix elements of the electron-impurity and e-e in- 
teractions, respectively, while cross terms are of higher 
order, e.g., Weel^'eiP- Thus, the two interactions give ad- 
ditive contributions to the conductivity. In Ref.^^, a phe- 
nomenological equation of motion was used to find the 
spin current in a system with disorder and longitudinal 
spin-Coulomb drag coefficient. We can use a similar ap- 
proach to obtain transverse spin currents with transverse 
spin-Coulomb drag coefficient l/r'j^- In the circularly- 
polarized basis. 



i{uj =F wo)jd 
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and correspondingly the spin-conductivities are 
n 1 



4m,., -(cj =F t^o)« + iM'" + l/rj;^ 
In the dc limit, this gives 

_ n 1/rf" + 1/rf 



a AO) 



4™* uj^+ {l/Tf"" + i/Ti^y 



(15) 



(16) 



Using Eq. , an identification of the e-e contribution is 
possible in a perturbative regime where l/rf^^ ^ 

luq, leading to the following formula: 
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(17) 



Comparison with Eq. (|13|1 enables us to immediately 
identify the microscopic expressions for the two scatter- 
ing rates. For the disorder contribution, we recover what 
we already knew, i.e., Eq. (jlip . For the e-e interaction 
contribution, we obtain 
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— lim ■ 

nm, V uj^Q 



,(18) 



where F'-^ is just the Coulomb force, and the force-force 
correlation function is evaluated in the absence of disor- 
der. The correlation function in Eq. ([TS]) is proportional 
to the function F^_{u>) which appeared in Ref. ^ [Eqs. 
(18) and (19)] in a direct calculation of the transverse 
spin susceptibility. Making use of the analytic result for 
'^mF^^ [ljj) presented in Eq. (21) of that paper we obtain 



Sao T^rtm^ajkl 
T'^l 27 (H-p)i/3 



(19) 
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FIG. 1; (Color online) The Gilbert damping q as a function 
of the disorder scattering rate 1/t. Red (solid) line shows the 
Gilbert damping for polarization p = 0.1 in the presence of the 
e-e and disorder scattering, while dashed line does not include 
the e-e scattering. Blue (dotted) and black (dash-dotted) lines 
show Gilbert damping for p = 0.5 and p = 0.99, respectively. 
Wetookg = O.lfcF, T = 54K, uo = EF[{l+pf^^~{l-pf^^], 

sm~^, rs — 5, a, = 2ao 



(14) Mo = 7pn/2, m* = rUe, n = 1.4 x 10^ 
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where T is the temperature, p = (n-f — ni)/n is the degree 
of spin polarization, a» is the effective Bohr radius, is 
the dimensionless Wigner-Seitz radius, ao = (4/97r)-'^/'^ 
and r(p) - a dimensionless function of the polarization 
p - is defined by Eq. (23) of Ref. [ill. This result is valid 
to second order in the Coulomb interaction. Collecting 
our results, we finally obtain a full expression for the 
Gilbert damping parameter: 

One of the salient features of Eq. (PO)) is that it scales 
as the total scattering rate in the weak disorder and 
e-e interactions limit, while it scales as the scattering 
time in the opposite limit. The approximate formula 
for the Gilbert damping in the more interesting weak- 
scattering/strong-ferromagnet regime is 

" 4TO,c^2^folrf- rr/' 
while in the opposite limit, i.e. for uq <C 1/t^'^, l/r'j^: 

4m* Mo Vrf" ^IV ' 

Our Eq. ([20]) agrees with the result of Singh and 
Tesanovio^ on the spin-wave linewidth as a function of 
the disorder strength and wq. However, Eq. (|^ also 
describes the influence of e-e correlations on the Gilbert 
damping. A comparison of the scattering rates originat- 
ing from disorder and e-e interactions shows that the lat- 
ter is important and can be comparable or even greater 
than the disorder contribution for high-mobility and/or 
low density 3D metallic samples. Fig. 1 shows the be- 
havior of the Gilbert damping as a function of the dis- 
order scattering rate. One can see that the e-e scatter- 
ing strongly enhances the Gilbert damping for small po- 
larizations/weak ferromagnets, see the red (solid) line. 
This stems from the fact that l/r^'" is proportional to 
l/r and independent of polarization for small polar- 
izations, while ^/t'^ is enhanced by a large prefactor 



T{p) = 2A/(1 - A2) + (l/2)ln[(l + A)/(l - A)], where 
A = (1 -p)^/^ /{1+pY/^. On the other hand, for strong 
polarizations (dotted and dash-dotted hnes in Fig. 1), the 
disorder dominates in a broad range of 1 /t and the inho- 
mogenous contribution to the Gilbert damping is rather 
small. Finally, we note that our calculation of the e-e in- 
teraction contribution to the Gilbert damping is valid un- 
der the assumption of hw <^ ksT (which is certainly the 
case if a; = 0). More generally, as follows from Eqs. (21) 
and (22) of Ref. [ll|, a finite frequency lu can be included 
through the replacement {2TrkBT)^ {iirkBT f + {hw f 
in Eq. p^ . Thus I/t*^ is proportional to the scattering 
rate of quasiparticles near the Fermi level, and our damp- 
ing constant in the clean limit becomes qualitatively sim- 
ilar to the damping parameter obtained by Mineev^ for 
Lo corresponding to the spin- wave resonance condition in 
some external magnetic field (which in practice is much 
smaller than the ferromagnetic exchange splitting cjq)- 

Summary - We have presented a unified theory of the 
Gilbert damping in itinerant electron ferromagnets at 
the order , including e-e interactions and disorder on 
equal footing. For the inhomogeneous dynamics (q' ^ 0), 
these processes add to a q = damping contribution 
that is governed by magnetic disorder and/or spin-orbit 
interactions. We have shown that the calculation of the 
Gilbert damping can be formulated in the language of the 
spin conductivity, which takes an intuitive Matthiessen 
form with the disorder and interaction contributions be- 
ing simply additive. It is still a common practice, e.g., in 
the micromagnetic calculations of spin-wave dispersions 
and linewidths, to use a Gilbert damping parameter in- 
dependent of q. However, such calculations are often at 
odds with experiments on the quantitative side, particu- 
larly where the linewidth is concernedi^ We suggest that 
the inclusion of the q^ damping (as well as the associ- 
ated magnetic noise) may help in reconciling theoretical 
calculations with experiments. 
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